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ABSTRACT 

In the pursuit of a general formulation for a modified gravitational theory at the non- 
relativistic level and as an alternative to the dark matter hypothesis, we construct 
a model valid over a wide variety of astrophysical scales. Through the inclusion of 
Milgrom's acceleration constant into a gravitational theory, we show that very general 
formulas can be constructed for the acceleration felt by a particle. Dimensional ana- 
lysis shows that this inclusion naturally leads to the appearance of a mass-length scale 
in gravity, breaking its scale invariance. A particular form of the modified gravitational 
force is constructed and tested for consistency with observations over a wide range of 
astrophysical environments, from solar system to extragalactic scales. We show that 
over any limited range of physical parameters, which define a specific class of astro- 
physical objects, the dispersion velocity of a system must be a power law of its mass 
and size. These powers appear linked together through a natural constraint relation 
of the theory. This yields a generalised gravitational equilibrium relation valid for all 
astrophysical systems. A general scheme for treating spherical symmetrical density 
distributions is presented, which in particular shows that the fundamental plane of 
elliptical galaxies, the Newtonian virial equilibrium, the Tully-Fisher and the Faber- 
Jackson relations, as well as the scalings observed in local dwarf spheroidal galaxies, 
are nothing but particular cases of that relation when applied to the appropriate mass- 
length scales. We discuss the implications of this approach for a modified theory of 
gravity and emphasise the advantages of working with the force, instead of altering 
Newton's second law of motion, in the formulation of a gravitational theory. 

Key words: gravitation - galaxies: kinematics and dynamics - galaxies: general. 



1 INTRODUCTION 

The dynamical mass to light ratios derived for spiral galax- 
ies are usually much greater than expected for their stellar 
components. This is often interpreted as indicating the grav- 
itational dominance of hypothetical dark matter. Alternat- 
ively, one could argue that the discrepancy between dynam- 
ical mass and baryonic mass is telling us that the Newto- 
nian law of gravity is not the one governing the dynamics. 
In particular , the Mo d ified N ewtonian Dynamics (MOND) 
proposed by iMilgronil (|l983bl ) has been proven to be suc- 
cessful in explai ning how gal axies rotate, without any dark 
matter (see e.g. ISanders fc M cGaugli 200^, for a review) . 

Recently, the range of astrophysical problems treated 
under the MOND hypothesis has increased significantly. 
Abundant recent publications on velocity dispersion meas- 
urements for stars in the local dwarf spheroidal galaxies, 
the extended and flat rotation curves of spiral galaxies, 
the large dispersion velocities of galaxies in clusters, the 
gravitational lensing due to massive clusters of galaxies. 



and even the cosmologically inferred matter content for the 
universe, have been successfully modelled under MOND. 
These, not as indirect evidence for the existence of a dom- 
inant dark matter component, but as direct evidence for 
the failure of the current Newtonian and general relativistic 
theories of gravity, in the large s cale or low acceleration 



regimes relevant for the above (sec Zander s fc Noordermeer 
2007; Nipoti ct al. 2007; Famacy ct al. 2 00"7l: iGentile et al 
200';i: iTiret et all boOTi: ISanchez-Salcedo et al.l I200S 



BekensteinI |2004 Icapozzieflo et al.l l2007l: ISoboutil boO? 



Mendoza fc Rosas- Guevaral |2007| . for recent examples) 



MOND has proved successful on many ast rophysica l situ- 
ations , though difficult on others ( see e.g. IMilgronil |2008| . 
l2009bl : lBekenstein|[200^ : IZhaol l2005l . for a good review on 
these points). 

The key feature of Milgrom's Modified Newtonian Dy- 
namics is the introduction of a fund amental acceler ation 
scale ao = 1.2 x 10"^°ms"^ (see e.g. iMilgromI [2008h into 
gravitation. The introduction of ao can alternatively be re- 
garded as grounded upon direct empirical evidence, as the 
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observed dynamics of large spiral discs attest. Additionally, 
non-relativistic gravity due to a point mass M, results in a 
force on any test particle such that it is pulled towards M 
with an acceleration a. Fundamentally, Newton's constant 
of gravity G completes the description of the problem and 
so, by m eans of Buck ingham's theorem of dimensional ana- 
lysis (cf. ISedovlll959l ), the absolute value of the attractive 
acceleration felt by the test particle located at a distance r 
from the point mass M is given by 



where 



with, 



a = ao /(x). 



1/2 



(1) 



(2) 



(3) 



The acceleration expressed in equation ([T]) converges to 
Newton's gravitational acceleration when the function 
f{x) = x'^ and to MOND's acceleration when f{x) = x. 
These two examples of functions / represent the gravita- 
tional approach to an extended non-relativistic theory of 
gravity. The main problem is how to find a function f{x) 
which, for the appropriate limits, converges to the Newto- 
nian and MONDian regimes. 



[ Hernandez et al.l (|2010l ) showed that the iBekensteinI 
(I2OO4F function f{x) = x + x^ serves quite well when ap- 
plied to dwarf spheroidal (dSph) galaxies and to the rota- 
tion curves at large radii of spiral galaxies, but is inconsistent 
with measured limits on departu res from Newtonia n grav - 
ity at solar system scales. Also, iFamaev fc BinnevI (|2005l ) 
showed that this particular form of the function f(x) does 
not work well when applied to our Galaxy. Despite the fact 
that this prescription has the corresponding limits as expec- 
ted {f{x) X when x and f{x) x^ when a; ^ 1), a 
more general function must be constructed. 

Note that pure dimensional analysis, with the introduc- 
tion of an acceleration scale ao, determines exactly the di- 
mensional form that the acceleration must have. In very gen- 
eral terms, it also shows that the introduction of ao means 
that gravity has a characteristic mass-length scale lyi which 
makes possible the construction of equation With all 
these, the acceleration turns out to be a function of the 
variable x only, which as we will show later, gives a robust 
way of working with an extended theory of gravity at the 
non-relativistic level. 

It is important to note that Milgrom has alr ea dy in- 
troduce d the length Im ( see e.g . 'MilgromI Il983bl lal. 1 19861 : 
[Milgro m fc Sanders! I2OO8I ; I Milgrom. ,200S . where it appears 
as a transition radius). In these studies it is shown 
that this mass-length scale serves as a transition point 
where the MONDian regim e passes to the Newtonian one. 
I Milgrom fc Sanders! (|2008l ) stressed the points that a mass 
distribution whose length is much greater than its associated 
mass-length Im is in the MONDian regime (since a; <^ 1) and 
a mass distribution whose length is much smaller than its 
mass-length scale is in the Newtonian regime (since a; ^ 1). 
The case a; = 1 can roughly be thought of as the point where 



the transition from the Newtonian to the MONDian regimes 
occurs. 

We now show that there is a connection between this 
approach and the one commonly used in the implementation 
of the MO ND theory. The physical forra of MOND is g iven 
by (see e.g. iBekenstein fc Milgrom!ll984l : lMilgrom!l200ll ') an 
Aquadratic Lagrangian (AQUAL) and so, its variation re- 
produces the equation of motion. For relevant symmetries in 
the problem, this approach gives the important result that 
the absolute value of the acceleration felt by a test particle 
in the presence of a point mass is given by 



a fi (a/ao) = |V(;/>n| 



GM 



(4) 



In this equation, the Newtonian scalar potential is repres- 
ented by and the interpolation function fi{a/ao) is such 
that /i(a/ao) = 1 in the Newtonian limit, which corresponds 
to a ^ ao and /x(a/ao) = a/ao in the MONDian regime, 
with a <C ao. Equating the acceleration in relation with 
that of equation ([l]), it follows that 



^(a/ao) 



(5) 



which implicitly shows that the MOND formalism can be 
equivalently expressed through the modification of the grav- 
itational force IT]). We note here that the MOND formula- 
tion (j4]) refers to a modification of the dynamical sector of 
the theory, whereas equation ([T]) is completely based on the 
modification of the gravitational force. Both are operation- 
ally equivalent formulations. The MOND formulation has 
always been tackled through dynamical modifications. How- 
ever, we show in this article that there are many advantages 
when choosing the modification in the gravitational sector. 
Therefore we only use the constant ao for consistency with 
the dynamical modifications. It is important to emphasise 
that in the gravitational modifications it is more natural 
to frame the problem in terms of the mass-length scale Im 
defined in equation ((3}. Furthermore, it is the use of dimen- 
sional analysis which tells us the very important fact that 
the dimensionless force f{x) in equation ([T]) only depends 
on the ratio Im/r. 

The article is organised as follows. Section [2] introduces 
a particular form of the function f{x). This is used in the 
subsequent sections for applications in different astrophys- 
ical environments, from solar system to galaxy cluster scales. 
Finally, in Section|3]we discuss the advantages of such a gen- 
eral function f{x). 



2 THE FORCE MODEL 

The dimensionless gravitational force f{x) in equation ([T} 
felt by a given test particle must be analytic, and as such it 
can be written as 



(6) 



We now show how to obtain a reasonable f{x) by simple 
exploration of the Newtonian and MONDian regimes. First 
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of all, notice that in the Newtonian and deep MOND limits, 
the function f{x) is such that f{x) — cn x'^ and f{x) = 
cmx respectively, with cn = cm = 1. We now focus on the 
Newtonian x 1 regime and explore an expansion about 
that limit of the form 



S)-(t)" + (t) + (t)° + (t " 



+ x-'' + ...) 



(7) 



The limit a: — s- oo gives the Newtonian acceleration, so tak- 
ing into account all the terms of the geometric series for 
X > 1 then. 



ao J \ X — 1 J X — 1 



(8) 



We now put special emphasis on the MONDian x <^ 1 
regime and explore the corresponding expansion, given by 



ao/ \ r / \ r J \ r 



= a; (l + x + + i'' + . . .) . 

(9) 

Note that the deep MOND regime is obtained in the limit 
X ^ 0, and so the geometric series of equation ((9)1 for x < 1 
gives 



This satisfies the Newtonian and MONDian acceleration 
limits for a; ^ oo, respectively. Note also that the case 
n = 1 with a minus sign is the same as two times the case 
n = with a plus sign, and both corresp ond to the Beken - 
stein ground state acceleration formula jBekensteinl 120041 '). 
This has prove d to be useful for the d ynamical modelling of 



dSph ga laxies (iHernandez et al.l2010l ), but not for our own 
Galaxy (jPamaev fc Binnevll200 5lV 



The acceleration function (|12p has no singularities, since 
according to I'Hopital's rule, a/ao — >■ (n + l)/n as a; — J- 1. 
In fact, to see this directly, notice that for the minus sign it 
follows from equation (|12|l that 



(1 - a;) (1 + a; + + a;^ + . . . + a:" 
(1 - a:)(l + a; + a;2 + . . . + x"-^) 



(1 + a: + a:^ + a:^ + 



(13) 



(1+a: + a:2 + . . . + 



For further applications we note that the right hand side 
of equation \12\ with a minus sign can be Taylor expanded 
as follows: 



I n + l n + 2 . 2n + l 2n + 2 , .\ 

= X + X —X +X —X +..., (14) 

for X < 1, and 



1 \ 2 1-n I 2-n l-2n , 2-2n , , ^ r-\ 

— =a:— a: + x — x + x +..., (15) 
ao' 



for X > 1. Choosing the positive sign we obtain: 



O- \ n+l , n + 2 I 2n+l 2n + 2 , fic\ 

— =a; — a; -|-a; -|-a; — x +..., (16) 



for a: < 1, and 



1 - X 



(10) 



Equation ((TJ can be thought of as the series for the 
negative powers of relation ((6)1 and equation Q as the one 
for the positive powers of the same relation. The interesting 
thing to note is that both of them can be analytically con- 
tinued for all values of x. For the limit cases, the minus sign 
on the denominator of both equations ((8} and (|10|l can be 
changed for a positive sign. 

Since we are interested in the complete analytic series 
let us propose a general acceleration formula given by the 
addition or substraction of equations ^ and (|10|l as follows: 



X ± a; 
1 ±1 



(11) 



Note that this last equation tends to the Newtonian acceler- 
ation regime when a: ^ oo and to the MONDian acceleration 
limit when a; — >■ 0. In fact, due to the symmetry of the nu- 
merator and denominator of equation pip , a more general 
relation can be postulated: 



l±a:" 



l±x" 



(12) 



a \ 2 , \-n 2-n l-2n , 2 — 2n . 

= a:+a: — x — x -fa; 4-..., (17) 

ao ' 



for a; > 1. This shows that the Newtonian and MONDian 
regimes are reached in the correct limit regardless of the 
value of n. 

The extreme limiting case of n — ^ oo corresponds to the 
function 



ao 



for < a: < 1 
for a- > 1 



(MONDian regime). 
(Newtonian regime). 



(18) 

This acceleration formula is of no use due to the discontinu- 
ity on the first derivative at a; = 1, but serves as a reference 
to understand that the real acceleration must smoothly pass 
fr om the Newtonian to the M ONDian regime. Also, as noted 
bv lMilgrom fc Sanders! (|2008l ). the point where x — 1 repres- 
ents approximately the transition from the Newtonian to the 
MONDian regimes. In the proposed model, this is strictly 
valid in the extreme case, with n — >■ oo. This point is relev- 
ant, as any function moving away from the value a = oq at 
a; = 1 seems to have a better chance at modelling a more real 
astrophysical situation, since it will smoothly transit from 
the MONDian regime to the Newtonian one. 

Figure [T] shows a plot of a/ao as a function of x for 
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Figure 1. The figure shows the acceleration function a in units of 
Milgrom's constant ao as a function of the parameter x. The thick 
dash-dot curve is the extreme limiting value n — >■ oo, i.e. a/ao = x 
for X 4; 1 and a/ao = x^ for x ^ 1. The curves above and below 
this extreme acceleration line represent values of n = 4, 3, 2, 1, 
for the minus and plus signs of equation II12I I respectively. The 
extreme limiting curve has a kink at a; = 1 and is of no physical 
interest due to the undefined derivative at that point. 



various values of n and choice of signs in equation (|12|) . 
From this figure it is seen that the curves with n ~ 3 — 4 are 
very close to the extreme limiting case, but preserve a good 
soft transition region between the MONDian and Newrtonian 
regimes. However, it is through fits with observations that 
an optimal real number n is to be calibrated. In fact, the fit 
with observations must also give us a way to decide between 
the plus or minus sign in equation (|12p . 

In order to fix the parameter n and the choice of sign in 
equation (|12|) . we construct different corresponding MOND 
interpolation functions fj,{a/ao) for this acceleration formula 
and compar e them with the best m odel for our Galaxy 
presented bv lFamaev fc BinnevI (|2005l '). To do that, we must 
substitute equations p2|> and ^ into ©, yielding 



fi{a/ao) = : 



1 ± 



(19) 



The function x{a/ao) that appears in equation (|19[) is ob- 
tained by solving numerically equation (|12|l for a fixed value 
of n. Figure |2] shows that the best fit t o the optimal /^(a/ao) 
obtained bv lFamaev fc BinnevI (|2005l ') for our Galaxy, is re- 
produced with the minus sign and with n = 3.13 « 3. The 
effective gravitational acceleration formula is hence chosen 
as 



1 + X + 



(20) 



Given observational errors, and uncertainties in the mass to 
light ratios and their radial variability in our galaxy, the con- 
strains of Figure [2] are subject to considerable uncertainties. 
As such, we intend to show mearly that functions of the pro- 
posed family are clearly consistent with available estimates 
of fx{a/ao), in this rather uncertain parameter range. In the 
following section we show this particular f{x) as consistent 
with the much more stringent constraints available and the 



Figure 2. The figure shows MOND's interpolation func- 
tion fi(a/ao) for our G alaxy (dotted curve) as inferred by 
iFamaev fc BinnevI ll2005f) . From bottom to top (not count- 
ing the dotted cu rve) the plots represent the foll owing mod- 
els: ( 1) iBekensteirJ ground state acceleration model ijBekensteinl 
|2004) . which corresponds to our equation II12I I with n = 1 and 
a minus sign. (2) The in terpolation formula p-{x) = x/(l + x) 
llFamaev fc BinnevI [ioosl. dashed). (3) n = 2 for a minus sign. 
(4) n = 3 for a minus sign. (5) The interpolation formula 
^(x) =x/ [1 + X^)^^^ l|Milgromlll983bl . dashed). (6) n = 3, plus 
sign. (7) n = 2, plus sign. (8) n = 1, plus sign. 



quasi-Newtonian scale of the solar system. Also, simple ex- 
planations for the observed structural relations of elliptical 
and dwarf spheroidal galaxies will be shown to appear nat- 
urally in Section [3| Note that for this particular case, the 
function /i(a/ao) has an analytic solution, since the func- 
tion x{a/ao) from equation H12|) with n = 3 is the root of a 
fourth order polynomial in x. However, due to its complic- 
ated form, we omit it here. It is interesting that for this value 
of n, the expansion about a; ^ 1 begins with the Newtonian 
term, and then skips the following two terms according to 
equation (|15[l . This guarantees that dynamics will remain 
extremely close to Newtonian for a large range of values of 
X > 1 and sheds light on the extended Newtonian character 
of the force of gravity. Similarly, for x <^ 1, the leading term 
of equation (|14p gives the deep MOND regime x with the 
following one being a;*. The absence of the x'^ and a;^ terms 
implies that physics close to the deep MOND regime will 
not present any strong variations for a considerable range of 
values in a; < 1. 

We also note that the acceleration function ([ij is such 
that Newton's theorems, i.e. the acceleration field at dis- 
tance r from the centre of a spherical system depends only 
on the total mass M(r) interior to r, while external shells 
result in no force, are valid for any analytic function ((Gjl 
which depends on the parameter x only. In order to see this, 
suppose that a cx a;'' with p an integer number. Assume that 
the test particle is placed at position r inside a spherically 
symmetric shell. If we trace a cone with solid angle Si} and 
vertex at the test particle, the shell is intersected at two 
opposite points ri and r2- The masses SMi and 5M2 con- 
tained within the solid angle Sil at these points keep the 
proportion &M2/&M\ = (r2/ri)^. This relation means that 
{SMi/rlY^'^ = (<5A^2/ri)^'/^ and so 5xP{r-i) = 5xP{r2). In 
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other words, the acceleration exerted by the outer shell at 
position r cancels out. Since we can do this for any integer p, 
it follows that any analytical function of a:: (cf. equation ((6])) 
guarantees Newton's theorems. This is an encouraging prop- 
erty of the acceleration ^ . 

We now extend equation 11} to a vector form. As noted 
above, the AQUAL formulation is satisfied by our model and 
as such its f ield equation (Poisson's gene ralisation formula) 
is given by (|Bekenstein fc Milgromlll984l 'l 



V • [^(a/ao) V4>] = 47rGp : 



(21) 



where the scalar potential <j) satisfies the condition a = 
— For systems with a high degree of symmetry, 
iBekenstein fc Milgroriil (| 19841 ) showed that 



do 



M (a/fflo) 



ao 



-V(/)N 



G_ 

ao 



(22) 

where we can now express the interpolation function ^ 
through expression ([5}, and p{r) represents the matter dens- 
ity at the coordinate point r. On the other hand, the vector 
form of equation ((T} must necessarily be of the form 



aof{x)ea, 



(23) 



where the unit vector Sa points in the direction of the accel- 
eration a. Substitution of equation (|23|) into (|22|) with help 
of relation ((5]) yields 



G_ 

ao 



■p{r')dV'. 



(24) 



This description generalises the previous relations in the 
sense that a point mass Af can be directly substituted for 
M{r) in all relevant equations. 

The magnitude of the vectorial equation (|24|l is the for- 
mula to calculate the variable x for a given mass distribu- 
tion density p{r). As an example, for a point mass M the 
density is given by p{r) = MS{r), where S represents Dirac's 
delta function. Consequently, the variable x is given by equa- 
tion ((2)1 as expected. For the case of a spherically symmet- 
ric distribution of matter, the density depends only on the 
radial coordinate r. This means that (|Binnev fc Tremaind 
I2OO8I ) x'ea = -GM{r)er/r', where M(r) is the mass con- 
tained within the radius r and Sr is a unit vector in the 
direction of the radial coordinate. With this, it follows that 
the variable x depends only on the mass contained within 
the radius r. Since the general acceleration function f{x) 
is analytic (cf. equation (O), then it is also clear from this 
point of view that Newton's theorems are valid for the class 
of models pre sented in t his ar ticle. 

Recently iMilgromI (|2010l ) has developed a quasi-linear 
formulation of MOND, which in particula r can by applied 
to ou r spherical symmetric case (see also IZhao fc FamaevI 
|2010D . We now show that this is equivalent to our for- 
mulation. In this theory, a general gravitational potential 
"1> := 0N + is proposed, where satisfies the equation 



where v{y) represents a new interpolation function which 
tends to y~^^^ far away from the strongest gravity regime. 
Our approach is equivalent to their results, since the con- 
nection between f{x) and v is given by 



fix) = x^ [u{x^) + 1] 



(26) 



Comparing this last result with equation ([5} it follows that 
p. = (u + and so, for i/ 3> 1, then p = u"^. For the 
specific value i^iy) ~ y^^^^ then 



(27) 



The above two equations show that the function / is spe- 
cifically a function of x^ only, i.e. it depends exclusively on 
the Newtonian force x^. 

In the following sections, we discuss several astrophys- 
ical systems where the proposed gravity law (|20|l can be 
tested over a wide range of values of x, and explore the 
various predictions which emerge. We perform comparisons 
mostly through expected scalings related to velocity disper- 
sions, masses and equilibrium radii, derived from very gen- 
eral dimensional arguments. We consider mostly spherical 
mass distributions M{r), making use of the validity of New- 
ton's theorems for spherical symmetry already discussed. 



2.1 Solar system consistency and an illustrative 
rotation curve of our Galaxy 

As a first test of the proposed force law we compare the vari- 
ations with respect to Newtonian acceleration which would 
be introduced at solar system scales, to the exquisitely meas- 
ured upper limits on this quantity. 

For solar system scales with distances of between 0.1 
and ~ 50 AU and taking M = 1 Mq, we get 10^ < a; < 10"* 
(cf. Figure[7| • Due to these large values, we note immediately 
that the deviations from Newtonian dynamics are negligible. 
Figure [3] shows that our model lies well within the observa- 
tional upper limits for the departures from Newtonian radial 
ac celeration for the planet s in our solar system, as reported 
bv lSereno fc Jetzed (|2006l ). In consequence, under the pro- 
posed force law, MONDian type dynamics have no relevance 
at solar system sc ales. In light of o ur result s, we argue in fa- 
vour o f studies by Anderson et al.l (2002) ; Turvshev fc TothI 
(|2009l ): iToth fc Turvshev! (|2009l )7 that regard the anomaly 
of the recession velocity of the Pioneer probes as an effect 
of the thermal radiation of the spacecrafts (see however, 
lBekensteinll2006l : iMilgromI [2009al . for interpretations of the 
anomaly as a signature of MOND). 

As a particular example, we now turn to the rotation 
velocity curve of our Galaxy. The circular rotation velocity 
V{r) associated to equation PU)) is given by 



Vir) 



aoh 



1/2 



(28) 



V ip - 



(25) 



For the case of our Galaxy, the mass of the bulge and the 
mass of the disc as a functi on of the radial distan ce are 
respectively given by (see e.g. lAUen fc Santillanlll99ll ) 
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Figure 3. The curve shows the fractional modifications Sa := 
(a — ONt) /iNt our model introduces when compared to purely 
Newtonian acceleration a^t- This is compared with the upper 
limits on t he deviations from rad ial Newtonian acceleration as 
reported bv lSereno fc JetzeJ l l2006ll . The diagram shows that the 
modifications introduced by our model are always within the ob- 
served upper limits. From left to right the crosses represent the 
calculated values Sa at the distance of Mercury, Venus, Earth, 
Mars, Jupiter, Saturn, Uranus, Neptune and Pluto respectively. 



m(r)\, , = mr^ (a^ 

^ ^ I bulge V 



, 2\-3/2 

+ r ) 



m(r) 



where m 



27rE()r, 



+ 1 e 



-r I 



1.4 X 10^" Mq, a = 0.387kpc, Eo 



(29) 
(30) 
7 X 



10* M0 kpc ^ and = 3.5 kpc. The effective mass at a 
distance r from the centre of our Galaxy is then given by 



M(r) = m(r-)| +m(r) 



(31) 



and so, by substitution of equation (|3ip into (|28p the the- 
oretical rotation curve of our Galaxy, shown in Figure IH 
is obtained. Given the validity of Newton's theorems for 
spherical mass distributions in our case, the inner rotation 
curve of Figure |4] is fully self consistent. This occurs since in 
that region, the dynamics are completely dominated by the 
galactic bulge, which has the same spherical symmetry as is 
assumed in the calculation. The same holds to a large de- 
gree for the calculated rotation curve beyond a few disc scale 
radii, where the mass distribution has essentially converged. 
In the intermediate regime, the assumption of spherical sym- 
metry in the calculation necessarily introduces some error. 
However, by analogy with the Newtonian case, where the 
rotation curve of an exponentially decreasing infinitely thin 
disc deviates only about 20% from tha t of the correspond- 
ing s pherical mass distribution (see e.g lBinnev fc Tremaind 
l2008l ). this error is probably small. The validation of this 
analogy through a full 3D implementation of our model, lies 
beyond the scope of this first presentation of our model. In 
the light of this, Figure [4] simply provides a graphical rep- 
resentation of the fact that our j(x) w as calibrated from 
our G alaxy's best fit iJ,{a/aQ) as given bv lFamaev fc BinnevI 
ll2005ll . 



Figure 4. The figure shows the rotation curve (continuous line) of 
our Galaxy using the theory of gravity proposed in this article (cf. 
equation | |20| |V We have assumed that only disc (dashed-dotted 
line) and bulge (dotted line) contributions to the rotation curve 
are present. 



This rotation curve represents well the observed 
features of our Galaxy's rotation curve (see e.g 
iKuiiken fc Gilmord Il989h . This was actually built into 
our model, as the value of n in the proposed force law 
was calibrated from the numerical MOND i nterpolation 
function estimated bv lFamaev fc BinnevI (|2005h . which was 
obtained precisely from calibrations against the rotation 
curve of our Galaxy. 



3 EQUILIBRIUM RADII OF BOUND 
CONFIGURATIONS 

If we think of an arbitrary astrophysical equilibrium struc- 
ture with total mass M, with characteristic equilibrium ra- 
dius and internal velocity dispersion a, we can estimate 
the relation between the above three quantities by equating 
the kinematic pressure to the mean gravitational force. Note 
that in this section we use M to denote the total mass of an 
extended astrophysical system. Writing the kinematic pres- 
sure per unit mass as Cpa^ /M, replacing p for A'lr^'^ and 
equating this to the proposed force law of equation (|20p we 
obtain 



1 



H,Q — Re 



(32) 



where we have introduced a dimensionless velocity disper- 
sion (T ~ a/{GaoMy'''^ and an equilibrium radius Rc := 
Tc/Im for the problem. The constant C is expected to be of 
order 1. 

We notice that for small values of _Rc, corresponding 
to large values of the variable x, i.e. the Newtonian regime, 
the right hand side of equation (|32|l reduces to 1/Rc- In 
this case equation (|32} leads directly to the virial equilib- 
rium relation of Newtonian gravity rc = C^'''^a /{Gp)^'''^ . At 
the opposite MONDian limit, with small values of x and so, 
large values of i?e , the right hand side of equation (|32|> tends 
to 1. In this limit, equation (|32p yields the Faber- Jackson re- 
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lation C^^^a — (GaoAf)^/*. Allowing for a proportionality 
between isotropic velocity dispersions in pressure suppor- 
ted systems and rotation velocities in angular momentum 
supported ones, this last relation could also be seen as the 
baryonic TuUy -Fisher relation wh ich also has an index close 
to 4 (see e.g. iPuech et al.]|2O10D . It is interesting that in 
terms of the equilibrium gravitational radii of the proposed 
force law, the galactic Faber- Jackson relation appears as the 
"Jeans mass" solution for large values of Re- Intermediate 
cases should be well described by the appropriate points 
along the transition between the two limits of equation (|32p . 

If we think of a gravitationally bound system with an 
associated dispersion velocity cr, the above results can also 
be easily understood from the point of view of dimensional 
analysis. Indeed, employing the same procedure we used in 
constructing equation {T]), the problem is now characterised 
by the dispersion velocity a of the system, the radius r, 
the total mass M, Newton's constant of gravity G and Mil- 
grom's constant ao- With these, Buckingham's theorem of 
dimensional analysis demands that the velocity dispersion 
must have the following form 

cr = g{x), (33) 

with g{x) an arbitrary function to be determined. A natural 
approximation often used, is to model observations over a 
limited parameter range through power law representations. 
As such, we can explore the consequences of imposing a 
power law for the function g{x), i.e., 

a oc x"' — r^™. (34) 

By comparing equations (|32p and (|34[) . we see that 
equation (|34|) is merely a power law approximation to the 
full relation supplied by the proposed force law. One is 
lead to expect that over any limited range of values of x, 
equation (|32|) should be accurately approximated by equa- 
tion (|34p. with a suitable choice of the parameter q, which 
will more generally be a function of x. 

The strong prediction of equation (|20p . as seen in equa- 
tion (|34p . is that over any constrained range of values of x, 
e.g. corresponding to any well defined class of astrophysical 
objects, a relation of the type 

a = Cr^M^, (35) 

will always appear. Comparison with equation H34[) yields 
directly /3 — —a and shows that the two power law indices 
of equation (|35|) are not independent, but will obey the ne- 
cessary constraint 



Again, we see that the two limits explored above cor- 
respond to a = 1/2 and a = 0, for the Newtonian and deep 
MOND regimes respectively. It is clear that the constraint of 
equation (|36|) . obtained through dimensional analysis, is sat- 
isfied in these two cases. Note that this constraint can also 
be derived through a power law approximation to the full 
force model of equation (|32|) . In what follows it will become 



evident that this power law approximation and the corres- 
ponding predicted constraint l(36)l accurately reproduce the 
empirical scalings observed in elliptical galaxies and the well 
studied dwarf spheroidals of the local group. We shall also 
present an exact expression for ce(x) at all scales as a strong 
testable prediction. 

3.1 The fundamental plane of elliptical galaxies 

The clearest correlations between the structural parameters 
of elliptical galaxies are expres sed by the fundament a l plan e 
relations, first described by iDiorgovski fc DavisI |l983); 
iDressler et al.l (| 19871 ). as scaling relations between the ob- 
served effective radius r of the Galaxy, the central line of 
sight velocity dispersion cr, and the m ean surface brightness 
I within the effective radius (see e.g. iBernardi et al.l l2003l : 
iDesroches et al.ll2007l . for recent observations). Assuming a 
constant mass to light ratio and using the fact that the mean 
surface brightness is given by the luminosity L within the 
effective radius, such that / = L/2r^, measurements over 
large samples, now in the thousands of galaxies, give a rela- 
tion identical to equation H35p . 

In what follows we present a treatment for elliptical 
galaxies in isolation. Many such systems reside in galaxy 
clusters where external field effects might have some relev- 
ance, e.g. in acc ounting for the escape of high velocity stars 
(|Wu et al.l I2OO8I ). However, the well established constancy 
of the basic scaling relations for ellipticals across different 
density environments, implies that these external field ef- 
fects are by far not the physical causes driving the structure 
of the fundamental plane. As such, we shall not consider 
such environmental effects further at this point. 

A broad agreement in the literature y i elds v alues close 
to the ones reported by IBernardi et"aLl (|2003l ') of /3 = 
-0.3356 ±0.023 and 7 = 0.503 ±0.025. The first interesting 
point to note is that the above measured values deviate for 
P, slightly but noticeably, from the Newtonian expectations 
of /3 = —1/2, 7 = 1/2. This feature is termed the tilt in the 
fundamental plane, and has traditionally been explained in 
terms of mass to light ratios which vary as a function of 
mass for elliptical galaxies. Although such explan ations are 
not unreasonable (see e.g. IDesroches et al.ll2007l ). we shall 
see that this tilt in the fundamental plane can be accounted 
for, at constant mass to light ratios, by the slight deviations 
from the Newtonian limit expected from equation H32p at 
the values of 7?c which correspond to elliptical galaxies. 

As already seen from equation p2p under the power 
law approximation of equation H35p . we can compare the 
observed values for the power law scalings of the funda- 
mental plane of equation (|35p against the theoretical ex- 
pectations of our model through equation H36p . For the ob- 
served value of /3 = —0.3356 ± 0.023, equation p6p implies 
7 = 0.418 ± 0.012. This is marginally consistent with recent 
measurements of 7 = 0.503 ± 0.025 at a two sigma level, un- 
der the hypothesis of a constant average mass to light ratio 
for elliptical galaxies. We therefore see that a power law ap- 
proximation through a dimensional analysis approach to the 
problem furnishes the relation between the indices /3 and 7 
in equation H35p . given by relation (|36p . The actual value of 
both indices will be obtained theoretically through the use 
of the full force model in what follows. It is noteworthy that 
practically all of the tilt in the fundamental plane can be 
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Figure 5. The plot shows the variation of the power index a with 
respect to a;, in the general gravitational equilibrium relation (see 
text) valid for all systems and at all scales. For large values of 
the parameter x, the Newtonian limit a = 1/2 is reached, and for 
small ones convergence to the MONDian one « = is obtained. 
As shown in Figure [7] many astrophysical systems of interest fall 
in the intermediate regions about x = 1. 



range of physical parameters, the power law approximation 
to equation p2p ceases to be valid, and a warp will necessar- 
ily develop. Being elliptical galaxies very close to the Newto- 
nian limit of Q = 1/2, we can estimate the first order trend 
of this warp against mass by considering that typical dens- 
ities for elliptical galaxies tend to drop as one goes to more 
massive systems, which would move the fit further away from 
the Newtonian regime, leading to a gradual decrease in a 
with increasing total mass. By comparing with Figure 8 of 
iGargiulo et~alT(|2009l ). we see precisely this very trend, with 
the measured value for their index a in the log = a log (jq 
fit (notice that a oc a~^) increasing as the limit mass of 
the galaxies included in the sample increases. Therefore, the 
measured first order trends for the fundamental plane in- 
dexes with mass, are seen to agree with the expectations of 
the model presented. 



3.2 Local dwarf spheroidal galaxies 

We now turn to the well studied dSph galaxies of the local 
group and their scaling relations. Taking this time the first 
order terms for 7?o 3> 1 in equation (|32p . we obtain 



accounted for directly by the expectations of equation (|36p . 
The fullest explanation might well include a certain degree 
of mass dependence for the average mass to light ratios, 
as po inted out by various studies (see e.g IPesroches et al.l 
I2OO7I ). 

In order to evaluate the parameter a we proceed by fit- 
ting a power law function of the form (|34[l to equation (|32p . 
This can be done locally to the smooth function (|32p by 
calculating the tangent line to the function logcr(Iog2;) at a 
given point Po{xq,(jq), which satisfies the relation 

log G = a(x-o) log X + log (70 — a{xo) log xo. (37) 
The slope of the tangent line is given by 



a(x) 



dloga _ X 
9 log; 



Z + 2x + x^ 



2 {1 + x + x'^){l + X + x"^ + x-^ 



(38) 



which simply means that the constant a depends on the 
scale xo of the system as shown by Figure [5] Notice that 
a — > 1/2 for the Newtonian 2; — ^ 00 regime and a — >■ for 
the MONDian x -> limit. 

For any given class of astrophysical objects, the para- 
meter Q can now be obtained by evaluating x and then solv- 
ing for a from equation (|38|l . 

Note that typical values for elliptical galaxies with radii 
of order a few kpc, and masses one order of magnitude 
around IO^^Mq (see Figure [Tjl, we get a; ~ 1, which from 



equation (|38p yields a ~ 0.3. The tilt in the fundamental 
plane is hence naturally explained by the force law presen- 
ted, requiring only minimum departures from constant av- 
erage values for the mass to light ratios of elliptical galaxies. 

Recent studies spanning larger ranges of velocity 
dispersion and ma s s values for elliptica l galaxies (e.g. 
iGargiulo et~aLll2009l : iDesroches et al.ll2007l 'l. have begun to 
identify slight trends in the power indexes of the funda- 
mental plane. This is exactly what one would expect from 
the developments presented here, as when one broadens the 



l-f 



(39) 



Again we see that the first deviations from the deep MOND 
regime occur after skipping the and terms. This 
implies that near the deep MOND regime, gravitational 
physics remain largely unaltered. If we now solve for rc from 
the above equation, we obtain 



(GMao) 



al \Cg^ - (GMao)i/2]- 



(40) 



Through full integration of isothermal equilibrium dens- 
ity profiles for the local dSph's under a force law equi- 
valent to_wlm^_we_^ropose here in the deep MOND re- 
gime, [Hernandez et all (|2010l ) showed that a suitable de- 
scription of the problem can be found without the need to 
invoke the presence of any dark matter. Still, this might 
require the consistent inclusion of the disruptive effects 
of g alactic tides on the local dwa r fs, as pointed out b y 
Hernandej (j2007l ): lAngusI (|2008h : 



e^ 



Sanchez- Salcedo . 



McGaugh fc Wolj (l2010f) . It is interesting to remark that 
Hernandez et al. (|2010h also showed that there is a very 



clear correlation between the mass to light ratio resulting 
from the dynamical modelling in MOND type prescriptions 
(or alternatively, the dark matter fraction for these systems), 
and the relative ages of the stellar populations present, as 
directly inferred from statistical stu dies of their resolved 
stellar populations (see e.g. [Hernandez et al.l 2000l: DolphinI 
l2002l : lHelmill2008l : iMartin etaLll2008l : iTolstov et all 120091 '). 
This correlation is natural in a scenario where only the stars 
present act as a source of gravity, shining less as they age, 
but a ppears contrived under the dark matter hypothesis. 
Also, iHernandez et all (|2010l ) showed the expected scaling 
on average between a* and total stellar mass for the ob- 
jects in question, within observational errors, expected from 
the deep MOND regime. Introducing this scaling into equa- 
tion (1401) we obtain 
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Figure 6. Measured values of half light radii ri/2 S'lid internal 
velocity dispersions a for the local dSph galaxies with their as- 
sociated one sigma uncertainties. The solid line gives a best 
fit (T^ oc ri/2 relation, the exponent of which res ults from the 
propo sed force law. Data from the compilation of IWalker et al.l 
1120101 ). The point for the heavily distorted Sagittarius dwarf 
shown at the right end has been excluded from the fit. 



Te oc cr . (41) 

We can now test the expectations of our proposed force 
law through equation (|4ip . by plotting the measured val- 
ues of the velocity dispersion for the local dSph galaxies, 
against their observed half-light radii. We assume that our 
Tc in equation pif) will be equal to some constant times the 
measured half light r adius ri/2 and take t he compilation for 
these values found in IWalker et all (|2010l ') to plot Figure O 
We see that in spite of the large errors present, the expec- 
ted scahng of oc ri/2 given by the hne shown, clearly 
holds rather well. We have preferred this particular plane to 
test our model, as velocity dispersions and half light radii 
are the directly measured observables, not affected by the 
assumption of rather uncertain mass to hght ratios. 

Finally, we can try a first order extrapolation from the 
smallest Rc <^ 1 systems, the dSph galaxies, to the largest 
astrophysical structures in that regime, galaxy clusters, by 
taking values of r = ikpc and a = lOkms^^ from Figure [S] 
and increasing the velocity dispersion by a factor of 50. The 
expectation of equation (|41|l then being radii of 2.5 x 10^ 
times larger. Indeed, 2.5 Mpc and 500kms~ are represent- 
ative values for the structural parameters of large clusters of 
galaxies. Explaining the dynamics of galaxy clusters in the 
pure M OND formalism requires some un seen dark matter 
(see e.g. ISanderdlioO^ : lAngus et al]|201(]| '). so we defer the 
details of this problem to a latter study. Here we only point 
out that the first order scaling of equation (|4ip . extrapol- 
ated across ~ 6 orders of magnitude in radius, is at least 
not significantly in error. 

Again, we see that having a well defined force law, 
rather than dealing with a cumbersome MOND interpol- 
ation function, facilitates the physics considerably. This al- 
lows a more direct tracing of the astrophysical consequences, 
and yields more transparent predictions. Notice that the 
general results of the whole Section [3] are not dependent on 




1 10^ lO-" 10" 10** 10^" 10^^ 10^' 
M(Mo) 

Figure 7. The figure shows different values of x oc 
for typical masses M and radii r of different astrophysical sys- 
tems: (1) Solar system at Earth's orbit. (2) Solar system at 
Jupiter's orbit. (3) Solar system at the Kuiper belt radius. 
(4) Globular clusters. (5) Elliptical galaxies and bulges of spirals. 
(6) Dwarf spheroidal galaxies. (7) Outer regions of spiral galaxies. 
(8) Galaxy clusters. 



the details of /(x) , and will hold for any such function within 
the theoretical framework introduced. We note however, the 
good agreement with observations for the particular f{x) 
used. 

As a final summary of this section, we present Figure [71 
with typical values of the parameter x for average masses 
and radii for a range of astrophysical objects. This Figure, 
when compared to Figure [5] clearly complements the results 
of the previous sections. We see that solar system dynam- 
ics will be utterly unaffected by the modification proposed. 
Marginal modifications will appear in the dynamics of glob- 
ular clusters (to withi n observational uncertainties, see e.g. 
[Hernandez et aLlbOld ). A slight tilt is to be expected in the 
fundamental plane of elliptical galaxies and we see that the 
external regions of spiral galaxies, galaxy clusters and local 
dwarf spheroidal galaxies should appear as the most heavily 
"dark matter dominated" systems. 
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4 CONCLUDING REMARKS 

We have shown in this article that if a modification of either 
the dynamical or the gravitational part of Newton's grav- 
ity is to be performed, a gravitational modification has a 
greater advantage for physical clarity and ease of calcula- 
tions. Furthermore, it has become clear that with the addi- 
tion of Milgrom's acceleration constant ao, all physical rela- 
tions must be at least functions of the parameter x defined 
in equation ((2]), according to Buckingham's theorem of di- 
mensional analysis. As examples, we have shown that the 
acceleration felt by a particle on a centrally symmetric grav- 
itational field satisfies this condition. Also, we show that the 
velocity dispersion associated to a system in the linear ap- 
proximation gives rise to a generalised gravitational equilib- 
rium relation which converges to the well known results of 
the fundamental plane in the description of elliptical galax- 
ies, the TuUy-Fisher relation for spiral galaxies, the Faber- 
Jackson relation for elliptical galaxies, the Jeans' stability 
criterion in the Newtonian limit and to corresponding ob- 
served relations for dSph's galaxies. 

The inclusion of ao in the theory of gravity means that a 
characteristic mass-length has been added. This scale ap- 
proximately indicates the transition between the MONDian 
and Newtonian regimes. By a suitable choice of model, this 
transition can be built smoothly. On very general grounds, a 
characteristic scale Im has been introduced into gravitation 
and so, gravity is no longer scale invariant. 

We note also another testable prediction: as precision 
increases at solar system measurements the upper limits to 
departures from Newtonian acceleration will converge to def- 
inite values, as implied by Figure [S] 

The particular acceleration function we built is perhaps 
only an approximation to a more general gravitational law. 
The functions we built, although logical and quite precise in 
principle, may just be good candidates for the real gravita- 
tional theory and they should be thought of as such: good 
approximations which fit a wide variety of available astro- 
physical data to within errors. However, it has to be through 
an extension of the general theory of relativity that a more 
precise and fundamental form is to be found. This exten- 
ded relativistic theory of gravity must show that the PPN 
approximations converge in some limit to the acceleration 
function built in this article. 

The Newtonian character of gravity has not been con- 
vincingly proven across the astrophysical scales explored in 
this article and as such, modifications to the theory are feas- 
ible. It is also important to note that a very wide variety of 
modified gravity theories of the type developed here, includ- 
ing MOND, are fundamentally falsifiable: a single a > ao 
system appearing as heavily dark matter dominated, would 
invalidate them all as an alternative to the dark matter hy- 
pothesis. 
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